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^ ■ Abstract 
> 

^ . We present an explicit formula for the discrete power function introduced by Bobenko, 

which is expressed in terms of the hypergeometric r functions for the sixth Painleve equation. 
The original definition of the discrete power function imposes strict conditions on the domain 
^ . and the value of the exponent. However, we show that one can extend the value of the exponent 

^ I to arbitrary complex numbers except even integers and the domain to a discrete analogue of 

the Riemann surface. Moreover, we show that the discrete power function is an immersion 
when the real part of the exponent is equal to one. 



^ : 1 Introduction 



The theory of discrete analytic functions has been developed in recent years based on the theory of 
circle packings or circle patterns, which was initiated by Thurston's idea of using circle packings 
as an approximation of the Riemann mapping [18]. So far many important properties have been 
established for discrete analytic functions, such as the discrete maximum principle and Schwarz's 
lemma [6], the discrete uniformization theorem [15], and so forth. For a comprehensive introduc- 
tion to the theory of discrete analytic functions, we refer to [17]. 

It is known that certain circle patterns with fixed regular combinatorics admit rich structure. 
For example, it has been pointed out that the circle patterns with square grid combinatorics intro- 
duced by Schramm [16] and the hexagonal circle patterns [5, 8, 9] are related to integrable systems. 
Some explicit examples of discrete analogues of analytic functions have been presented which are 
associated with Schramm's patterns: exp(z), erf(z). Airy function [16], z^, log(z) [4]. Also, discrete 
analogues of and log(z) associated with hexagonal circle patterns are discussed in [5, 8, 9]. 
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Among those examples, it is remarkable that the discrete analogue of the power function 
associated with the circle patterns of Schramm type has a close relationship with the sixth Painleve 
equation (Pvi)[7]. It is desirable to construct a representation formula for the discrete power 
function in terms of the Painleve transcendents as was mentioned in [7]. The discrete power 
function can be formulated as a solution to a system of difference equations on the square lattice 
(n, m) e 7} with a certain initial condition. A correspondence between the dependent variable of 
this system and the Painleve transcendents can be found in [14], but the formula seems somewhat 
indirect. Agafonov has constructed a formula for the radii of circles of the associated circle pattern 
at some special points on 1} in terms of the Gauss hypergeometric function [3]. In this paper, we 
aim to establish an explicit representation formula of the discrete power function itself in terms 
of the hypergeometric t function of Pyi which is valid on = {(n, in) e 'I?\n,m > 0} and for 
7 6 C\2Z. Based on this formula, we generalize the domain of the discrete power function to a 
discrete analogue of the Riemann surface. 

On the other hand, the fact that the discrete power function is related to Pyi has been used 
to establish the immersion property [ '] and embeddedness [ ] of the discrete power function with 
real exponent. Although we cannot expect such properties and thus the correspondence to a certain 
circle pattern for general complex exponent, we have found a special case of Re y = 1 where the 
discrete power function is an immersion. Another purpose of this paper is to prove the immersion 
property of this case. 

This paper is organized as follows. In section 2, we give a brief review of the definition of the 
discrete power function and its relation to Pyj. The explicit formula for the discrete power function 
is given in section 3. We discuss the extension of the domain of the discrete power function in 
section 4. In section 5, we show that the discrete power function for Re y = 1 is an immersion. 
Section 6 is devoted to concluding remarks. 



2 Discrete power function 

2.1 Definition of tlie discrete power function 

For maps, a discrete analogue of conformality has been proposed by Bobenko and Pinkall in the 
framework of discrete differential geometry [10]. 

Definition 2.1 A map f : ^ C ; (n, m) is called discrete conformal if the cross-ratio 

with respect to every elementary quadrilateral is equal to -I: 

^fn,m fn+l,m)(^fn+l,in+l fn,m-^\) _ ^ ^2 
(,fn+l,m fn+\,m+l^i.fn,m+\ fn,m) 

The condition (2.1) is a discrete analogue of the Cauchy-Riemann relation. Actually, a smooth 
map / : Z) c C ^ C is conformal if and only if it satisfies 

lim ^-^^^^ ~ -^^^ ^ ^' + e,y + e)- f(x, y + e)) ^ ^ 
'^0 (fix + e,y)- fix + e,y + e)){f{x, y + e) - f(x, y)) 

for all (x,y) e D. However, using Definition 2.1 alone, one cannot exclude maps whose behaviour 
is far from that of usual holomorphic maps. Because of this, an additional condition for a discrete 
conformal map has been considered [2, 4, 7, 11]. 
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Definition 2.2 A discrete conformal map fn^m is called embedded if inner parts of different ele- 
mentary quadrilaterals ifn,m, fn+\,m, fn+i,m+i, fn,m+i) do not intersect. 

An example of an embedded map is presented in Figure 1 . This condition seems to require 
that / = fn m is a univalent function in the continuous limit, and is too strict to capture a wide 
class of discrete holomorphic functions. In fact, a relaxed requirement has been considered as 
follows [2, 4]. 

Definition 2.3 A discrete conformal map f^^m is called immersed, or an immersion, if inner parts 
of adjacent elementary quadrilaterals (fn^m, fn+i,m, fn+i,m+i, fn,m+i) <^re disjoint. 

See Figure 2 for an example of an immersed map. 

Let us give the definition of the discrete power function proposed by Bobenko [4, 7, 1 1]. 

Definition 2.4 Let f : ^ C ; (n, m) i-^ /„ „, be a discrete conformal map. If fn,m is the solution 
to the difference equation 

r ^ ifn+l,m ~ fi,ni)ifn,m ~ fn-l,m) ^ ifn,m+l ~ fn,m){fn,m ~ fn,m-\) ,^ 

yfn,m = 2n — + 2m — (2.3) 

Jn+l,m Jn-\,m Jn,m+\ Jn,m—\ 

with the initial conditions 

/o,o = 0, /i,o=l, fo^=ey'''l^ (2.4) 
for < 7 < 2, then we call f a discrete power function. 

The difl'erence equation (2.3) is a discrete analogue of the diff'erential equation yf = z — for 

oz 

the power function /(z) = zJ, which means that the parameter y corresponds to the exponent of the 
discrete power function. 

It is easy to get the explicit formula of the discrete power function for m = (or n = 0). When 
m = 0, (2.3) is reduced to a three-term recurrence relation. Solving it with the initial condition 

/o,o = 0, /i,o = 1, we have 



fnfl 



21 + y \ \ 2k -y ' 



, (2.5) 

nl^ (n=2/.i), 



for n e Z+. When m = 1 (or n = 1), Agafonov has shown that the discrete power function can be 
expressed in terms of the hypergeometric function [3]. One of the aims of this paper is to give an 
explicit formula for the discrete power function yj, ,„ for arbitrary {n, m) e Z^. 

In Definition 2.4, the domain of the discrete power function is restricted to the "first quadrant" 
Z^, and the exponent y to the interval < y < 2. Under this condition, it has been shown that the 
discrete power function is embedded [2]. For our purpose, we do not have to persist with such a 
restriction. In fact, the explicit formula we will give is applicable to the case y e C\2Z. Regarding 
the domain, one can extend it to a discrete analogue of the Riemann surface. 



3 



Figure 1: An example of the embedded 
discrete conformal map. 



Figure 2: An example of the discrete con- 
formal map that is not embedded but im- 
mersed. 



2.2 Relationship to Pyi 

In order to construct an explicit formula for the discrete power function we will move to 
a more general setting. The cross-ratio condition (2.1) can be regarded as a special case of the 
discrete Schwarzian KdV equation 

^fn,m fn+\,m)^fn+l,m+\ fn.m+Y) Pn ,^ 

= — , (2.6) 

^fn+\,m fn+l,in+l^(^fn,m+l fn,m) Qm 

where p„ and q,„ are arbitrary functions in the indicated variables. Some of the authors have 
constructed various special solutions to the above equation [12]. In particular, they have shown 
that an autonomous case 

^fn,m fn+l,in)(^fn+\,m+l fn,m+Y) _ 1 ^2 y-j 

^fn+\,m fn+l,m+l^(.fn,m+l fn.iii) ^ 

where t is independent of n and m, can be regarded as a part of the Backlund transformations of 
Pvi, and given special solutions to (2.7) in terms of the r functions of Pyi- 

We here give a brief account of the derivation of Pyi according to [14]. The derivation is 
achieved by imposing a certain similarity condition on the discrete Schwarzian KdV equation (2.7) 
and the difference equation (2.3) simultaneously. The discrete Schwarzian KdV equation (2.7) is 
automatically satisfied if there exists a function v„ ,„ satisfying 

fn,m fn+l,m ~ ^ ^n,m^n+\,nf> fn,m fn,m+\ ^n,m^n,m+l- (2-8) 

By eliminating the variable we get for v„_,„ the following equation 

which is equivalent to the lattice modified KdV equation. It can be shown that the difference 
equation (2.3) is reduced to 

^ V„,um-Vn-I,,n ^ ^^ V„,^,i - V„,„,-i ^ ^ _ ^_^yn,n^ (2.10) 

with 7=1-1- Ijx, where /I 6 C is an integration constant. In the following we take A = iiso that 
(2.10) is consistent when n = m = Q and Vi_o + v_i,o i^O Vqj + Vq-i. 
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Assume that the dependence of the variable v„,„, = Vn,m{t) on the deformation parameter t is 
given by 

^ d Vn+\,m ^n-\,m i i \ 

- It— log V„,m = n ■ +Xn+m, (2. 1 1) 

where ;^^„+„, = Xn+m{t) is an arbitrary function satisfying ;^f„+,„+2 = Xn+m- Then we have the following 
Proposition. 

Proposition 2.5 Let q = ^„ ,„ = qn,m{t) be the function defined by q^^m = Then q satisfies 



VI 



d^q 1/1 1 1 \ / dq\^ /I 1 I \dq 



dt^ 2\q^ q-1^ q-tj\dtj \t ^ t - I ^ q - t) dt 



q(q- l)iq-t) 
2t\t- 1)2 



-<-+<- —+(1-0^) 



V2 1(^-1)2 ^ \q-tf 



(2.12) 



Kl, = ^il^-v + m- nf, Kl = h4i-v-m + nf, 
= ^{pi + v-m-n- 1)2, 6*^ = ^Qj. + v + m + n + 1)^, 



(2.13) 



where we denote v = (-1)™"^"//. 



Li general, Pvi contains four complex parameters denoted by /Coo, kq, ki and 9. Since n, m e Z+, 
a special case of Pyi appears in the above proposition, which corresponds to the case where Pyi 
admits special solutions expressible in terms of the hypergeometric function. In fact, the special 
solutions to Pvi of hypergeometric type are given as follows: 

Proposition 2.6 [ ] Define the function T^'ia, b, c; t) (c ^ Z, n' e Z+) by 

^^,(^^^^^.^^^1 det(^(a + /-l,Z, + ;-l,c;0)i,,,,. (n^ > 0), ^^.14) 



with 



(p(a, b, c; = Co F{a, b, c; 

r(c) 



r(a-c+ l)rib-c + 1) , 
+ci^ f(i^) '^(« - c + 1, Z7 - c + 1, 2 - c; 0. 



(2.15) 



Here, F(a, b, c; t) is the Gauss hypergeometric function, Y{x) is the Gamma function, and cq and c\ 
are arbitrary constants. Then 

^0,-1,0^-1,-1,-1 

^ = ,-1.-1.-^-1,0 (2.16) 

~n' ',i'+l 

with T^j/'"' = T„'(a + k + \,b + I + 2,c + m + \ \t) gives a family of hypergeometric solutions to Pyi 
with the parameters 

Koo=a + n', KQ = b-c+l+n', ki = c - a, 6 = -b. (2.17) 



We call T„/(a, b, c; t) or rf;/ '" the hypergeometric t function of P 



VI- 
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3 Explicit formulae 

3.1 Explicit formulae for and v„ 

We present the solution to the simuhaneous system of the discrete Schwarzian KdV equation (2.7) 
and the difference equation (2.3) under the initial conditions 

/o,o = 0, /i,o = Co, /o,i = cxf, (3.1) 

where y = 2r, and cq and ci are arbitrary constants. We set co = ci = 1 and t = e"'(= -1) to obtain 
the explicit formula for the original discrete power function. Note that T„'{b, a, c; t) = r^'ia, b, c; t) 
by the definition. Moreover, we interpret F{k, b, c; t) for k e Z>o as F(k, b, c;t) = and T{-k) for 
(-1)^ 

k e Z>o as T(-k) = -— ^■ 
k\ 

Theorem 3.1 For (n, m) e Z^, the function /„ „, = f„^mit) can be expressed as follows. 

( 1 ) Case where n < m (or n' = n). When n + m is even, we have 

= c.r'W^^^^^-" -ni-N,-r-N^l,-r,t) 
^' {-r+\)NTn{-N+\,-r-N + l,-r + 2;t) 

n + m 

where N = — - — and (u)j = u(u + 1 ) • • • (m + j - \) is the Pochhammer symbol. When n + m 
is odd, we have 

r _ r-n + ^)n-\ r^j-N + 1 , -r - N + 1 , -r, t) 

(-r+lh_,T,A-N + 2,-r-N + 2,-r + 2-ty ^^-"^ 

n + m + 1 

where N = . 

2 

(2) Case where n > m (or n' = m). When n + mis even, we have 

f - ,, (^ + Ik-i ^m{-N + 2, -r - + 1, -r + 2; 
(-r + l);v rmi-N +l,-r-N + 2,-r + 2;t) 

n + m 

where N = — — . When n + mis odd, we have 

(r + l)^^i T^i-N + 2,-r-N +l,-r+l;t) 
(-r + l)^_i T,n{-N +\,-r-N + 2,-r+\;t) 

n + m + 1 

where N = . 

2 

Proposition 3.2 For {n, m) 6 Z^, the function Vn,m = v,j,,„(0 can be expressed as follows. 

( 1 ) Case where n < m (or n' = n). When n + mis even, we have 

^ n ir)N T„i-N +l,-r-N +l,-r+Ut) 

~ \-r+ 1)n t„(-N +l,-r-N + 2,-r + 2;ty ' ^ 
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n + m 

where N = — - — . When n + mis odd, we have 
2 



., tA-N +\,-r-N + l,-r+\-t) 
'T„(-N + 2,-r-N + 2,-r + 2-ty 



n + m + I 

where N = . 

2 

(2) Case where n > m (or n' = m). When n + m is even, we have 



^^-^ {r)N T,n{-N +l,-r-N +l,-r+l;t) 



n + m 

where N = . When n + mis odd, we have 

2 



T,„i-N + 2,-r-N + 2,-r + 2; t) 

Vn,m = -Cot 2 — — , (3.9) 

T„,i-N + \,-r - N + 2,-r + l;t) 

n + m + \ 

where N = . 

2 

Note that these expressions are applicable to the case where r e C\Z. A typical example of 
the discrete power function and its continuous counterpart are illustrated in Figure 3 and Figure 
4, respectively. Figure 5 shows an example of the case suggesting multivalency of the map. The 
proof of the above theorem and proposition is given in the next subsection. 

Remark 3.3 Agafonov has shown that the generalized discrete power function /„ „„ under the 
setting ofco = Ci = \,t = e^'" (0 < or < tt) and Q < r < \, is embedded [3]. 

Remark 3.4 As we mention above, some special solutions to (2. 7) in terms of the r functions o/Pyi 
have been presented [12]. It is easy to show that these solutions also satisfy a difference equation 
which is a deformation of (23) in the sense that the coefficients n and m of (23) are replaced 
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Figure 5: The discrete power function 
witli 7 = 0.25 + 3.35/. 



by arbitrary complex numbers. For instance, a class of solutions presented in Theorem 6 of [H] 
satisfies 

(ao + a2 + a4)fn,m 

, ^^fn+l,m ~ fn,ni)ifn,m ~ fn-\,m) , ^(.fn,m+l ~ fn,iii)ifi,m ~ fi,m-l) 

= (n- Qfi) 7 Z~7 +a2 + a4- m) — , 

Jn+l,m Jn-l,m Jn,m+l Jn,m-l 

(3.10) 

where or,- are parameters ofPy\ introduced in Appendix A. Setting the parameters as (ao , i , 0^2 , 0^3 , «4) 
(r, 0, 0, -r + 1 , 0), we see that the above equation is reduced to (2.3) and that the solutions are given 
by the hypergeometric t functions under the initial conditions (3.1). 



3.2 Proof of the results 

In this subsection, we give the proof of Theorem 3.1 and Proposition 3.2. One can easily verify 
that /„„, satisfies the initial condition (3.1) by noticing To(a,b,c; t) = I. We then show that 
and Vn,m given in Theorem 3.1 and Proposition 3.2 satisfy the relation (2.8), the difference equation 
(2.3), the compatibility condition (2.9) and the similarity condition (2.1 1) by means of the various 
bilinear relations for the hypergeometric r function. Note in advance that we use the bilinear 
relations by specializing the parameters a, b and c as 

n + m 

a = -N, b = -r-N, c = -r+l, N = (3.11) 

when n + mis even, or 

n + m + 1 ^ , ^ 

a = -r-N+l, b = -N, c = -r+l, N= , (3.12) 

when n + mis odd. 

We first verify the relation (2.8). Note that we have the following bilinear relations 



(c - i)rr>'- =(c-b- i)tr:-:v'~''~' + bj^'xij^'-^ 
ic - i)r;''"''- vr^-i = (c-b- i)r°'->'V>'->'-^ + z^rr T-^-^■-^ 



(3.13) 
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/ _ 0,-1,-1 0,-1,0 _ -1,-1,-1 1,-1,0 _ 7 0,0,0 0,-2,-1 



(Z. - a + DrrV'-'-' = - c + l)T°r^'"r-''°'-> + (c - a)r°r>'- 



(3.15) 



for the hypergeometric r functions, the derivation of which is discussed in Appendix A. Let us 
consider the case where n' = n. When n + mis even, the relation (2.8) is reduced to 

-rT'^'''\'y^'' = Nr^y^'V!'"''' - (r + A^)r^^'^Vr^ ^'•'''^ 

where we denote 

^'"''"^ = T,A-N + h, -r-N + h, -r + t), (3.17) 

for simplicity. We see that the relations (3.16) can be obtained from (3.13) with the parameters 
specialized as (3.1 1). In fact, the hypergeometric r functions can be rewritten as 

T°'-i'-i =T„(a + 1,Z7+ l,c) = T„(-A^+ l,-r-A^+ l,-r+ 1) = Tf/'''", (3.18) 

for instance. When n + m is odd, (2.8) yields 



-r^^^i^^ = (-r + A^)rri^t^VI^^^^<'] - A^rfP^^lr^l'"] 



(3.19) 



which is also obtained from (3.13) by specializing the parameters as (3.12). Note that the hyper- 
geometric r functions can be rewritten as 



Tn~^''^ = Tn(.a + \,b+\,c) = Tn(.-r - N + 2, -N + \ , -r + \) 

= T„i-N +l,-r-N + 2,-r+l) = t|,' 



[1,2,1] (3-20) 



this time. In the case where n' = m, one can similarly verify the relation (2.8) by using the bilinear 
relations (3.14) and (3.15). 

Next, we prove that (2.3) is satisfied, which is rewritten by using (2.8) as 

r --L 

tnm fit - in 

-/-^ = — — + — —■ (3.21) 

n,m tt+l,m n~i,m n,m+l n,m~\ 



(3.22) 



We use the bilinear relations 

ir n' ^ ' n' + l fi'-l ;i +1 fi'-l' 

/ ,7 , / , IX 0,0,0 0,-1,-1 -1,-1,-1 1,0,0 , /7 , IN 0,-1,0 0,0,-1 

(a + b-c + n +1)t^;' T^; = ar^, r ; ' +(Z?-c+1)t; ' t„; ' , 
and 

0,0,0 -1,-1,-2 _ _.-l -1,-1,-1 0,0,-1 , -1,-1,-1 0,0,-1 

' H ' « n+1 n-1 " " ' 

0,0,0 1,-1,0 _ 0,-1,0 1,0,0 _ 0,-1,0 1,0,0 /o --(ox 

' m ' m ~ 'm ' m m+l m-1 ' \D.Z.J) 

^0,-1,-1-1,0,-1 _ _ -1,-1,-1 0,0,-1 , -1,-1,-1 0,0,-1 
' m <m — 'm+l ' m-l ' ' 



for the proof. Their derivation is also shown in Appendix A. Let us consider the case where n' = n. 
When n + m is even, we have 

" " n+L n—l n+l n—l ' /o 

^_[l,2,2] [1,1,1] _ »r_[0,l,l] [2,2,2] ^ [1,1,2] [1,2,1] y^'^"*) 

from the bilinear relations (3.22) by specializing the parameters fl,Z? and c as given in (3.1 1). These 
lead us to 

, „ _[1,2,2] [1,1,1] 

-1 + -1 ^c-h-'^'T^-- - 

"^n+l,™ ^ "^n-Um ^1 [0,1,1] [1,2,1] ' 

n+l n-\ /o 

[1,2,2] [1,1,1] (.o.^j; 

-1 , -1 -1 ^-f-+2 ^« 

V + V = —C t 2 

'^n,m+l '^«,m-l "-1 ^ [0,1,1] [1,2,1]' 

By using 

'n '« ~ ' 'n+l 'n-1 ^ 'n '« 

which is obtained from the first relation in (3.23), one can verify (3.21). When n + m is odd, we 
have the bilinear relations 



_[ 1,2,2] [0,1,0] _ _^-l_[0,l,l] [1,2,1] _[0,1,1]_[1,2,1] /o 

',1 <„ — ' '„+] ',,„] ',1 5 v-^-^'J-' 



.„^[2,2,2]^[i,2,i] ^ ^ A^)r:it'V[f:f ] + (r + - Dr^^^^:^, 

,^[2,2,2]^[1,2,1] = + ^ _ 1)^[U.1]^[2,3,2] ^ ^ ^)^n,2,2]^[2,2,l] ^ 



(3.27) 



from (3.22) with (3.12), and 

_[2,2,2] [1,1,0] _ _^-l_[l,l,l]_[2,2,l] [1,1,1] [2,2,1] ^ox 
'ji 'n ~ ' ',,+1 ',,-1 'n '« ' V^-^o; 

from the first relation in (3.23). These lead us to (3.21). We next consider the case where n' = m. 
When n + mis even, we get the bilinear relations 

_ [1,2,2] [1,1,1] ^ ^ [1,1,2] [1,2,1] ^1 [0,1,1] [2,2,2] 

■'^'m+1 'm-1 'm+1 'm-1 ' H 90'> 

„_[1,2,2] [1,1,1] _ AT [0,1,1] [2,2,2] AT_[1, 1,2] [1,2,1] (.J-^^; 

and 

[1,2,2] [2,1,2] _ [1,1,2] [2,2,2] _ [1,1,2] [2,2,2] orvx 

'm 'm ~ 'm 'm 'm+1 'm-1 ' ^J.JU; 

from (3.22) and the second relation in (3.23), respectively. By using these relations, one can show 
(3.21) in a similar way to the case where n' = n. When n + mis odd, we use the bilinear relations 



and 

_[1,2,1] [2,1,1] _ __[1,1,1]_[2,2,1] _[1,1,1]_[2,2,1] ooN 
'm 'm ~ 'm+1 'm-1 m m ' V-^-J^^ 

which are obtained from (3.22) and the third relation in (3.23), respectively, to show (3.21). 

We next give the verification of the compatibility condition (2.9) by using the bilinear relations 

/ N 0,-1,-1 -1,-1,0 7 0,0,0 -1,-2,-1 IN -1,-1,-1 0,-1,0 

(c - aX; T^,;^ - br^; t,,;^ = {t - 1)t„, t^;^{ , 

/ N. 0,-1,-1 -1,-1,0 , 0,0,0 -1,-2,-1 ^. IN 0,-1,0 -1,-1,-1 ^ ■ ^ 
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The derivation of these is discussed in Appendix A. We first consider the case where n' = n. When 
n + mis even, we get 



from the bilinear relations (3.33). Then we have 

( . [1.1,2] [0,1,1] 

I Vn,m + Vn+l,m+l - t (.T [i_2,2] [0,1,2]' 

(-r + 1)n+i t\' 

. . [0,1,1] [1,1,2] 

> _ 1^ UJW ',,+1 

v„,M + r v„+i,„,+i - r -(,r i j [1,2,2] 10,1,2] ' 

from which we arrive at the compatibility condition (2.9). When n + m is odd, we have 

^ [1,2,1] [1,1,2] ^ [2,2,2] [0,1,1] , J. [1,1,1] [1,2,2] 
■'^'« 'n+1 ^■'^'n 'n+1 ^^'n '„+i j 

^ [1,2,1] [1,1,2] ^ [2,2,2] [0,1,1] ^ . ^ [1,2,2] [1,1,1] 
-'^"'i 'n+1 ^''^''1 'n+1 V* J^^'n 5 



(3.34) 



(3.35) 



(3.36) 



from (3.33). Calculating t~-Vn,m + Vn+\,m+\ and Vn,m + v„+i^+i by means of these relations, we see 
that we have (2.9). In the case where n' = m, one can verify the compatibility condition (2.9) in a 
similar manner. 

Let us finally verify the similarity condition (2.1 1), which can be written as 

n I d nvn+i.,„ 

X - oA'»+m - t-r log Vn,,n = ■ (3.37) 

^ ^ Vn+l,m + Vn-l,m 

Here, we take the factor ;^„+,„ as Xn+m = r[(-\)"'^'" - 1]. The relevant bilinear relations for the 
hypergeometric r function are 

/-n _i_ .^wO.CO -..0,-1,-1 ^,.-u-i,-i,-Ui,o,o 
(D + n)T„ ■ T„ = at T^^y T^_j , 

m , 1 , 1\ 0,-1,-1 0,0,0 ru , i\ 0,-1,0 0,0,-1 /o oo\ 

{D + b-c+\)T,n -T,- =(Z?-c+1)t,^ • t„; • , (3.38) 
/n _i_ ^ _L ,^^...o,o,o 0,-1,-1 ^ -1,-1,-1.^1,0,0 

The derivation of these is obtained in Appendix A. We first consider the case where n' = n. When 
n + m is even, it is easy to see that we have 

J0,1,1] [2,2,2] 

n =-Nr' 'llli^Tiir (3.39) 

Vn+l,m + Vn-l,m tJ," V}, 

from the bilinear relation (3.24). We get 

{D + n)r^^'^l . r^>'" = -NrW^^^'.'^^^fK (3.40) 

from the first relation in (3.38) with (3.11). From this we can obtain the similarity condition (3.37) 
as follows. When n + m is odd, we have 

(D + n)r[f'2'2l . ri'^^^i] = -r\r + - 1)t[';Y\S"\ (3-41) 
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from the first relation in (3.38). This relation together with the first relation in (3.27) leads us to 
(3.37). Next, we discuss the case where n' = m. When n + mis, even, we have 

(D + A^)Ti!'^-^l.ri;^'^'UA^ri|'>'2lTi^2^^ (3.42) 

from the second relation in (3.38). Then we arrive at (3.37) by virtue of the second relation in 
(3.29). When n + m is odd, we get 

(D + r+ '^)rl'/-'^ ■ = (r + N - DT^t'^'r^^^, (3.43) 

from the third relation in (3.38). Then we derive the similarity condition (3.37) by using the second 
relation in (3.31). This completes the proof of Theorem 3.1 and Proposition 3.2. 



4 Extension of the domain 

First, we extend the domain of the discrete power function to Z^. To determine the values of fn^m 
in the second, third and fourth quadrants, we have to give the values of /_i,o and /o -i as the initial 
conditions. Set the initial conditions as 

/-l,0 = C2?^ /o,-l=C3^^ (4.1) 

where C2 and cs are arbitrary constants. This is natural because these conditions reduce to 

/i,o = l, hi=e''"-, Uo = e'''\ /o,-i=e^"'" (4.2) 

at the original setting. Due to the symmetry of equations (2.7) and (2.3), we immediately obtain 
the explicit formula of yj, ,„ in the second and third quadrant. 

Corollary 4.1 Under the initial conditions /o,i = cif and (4.1), we have 

f-n,m — fn,m\ CQi-^cifi'' ■> f-n-m — fn,m\c(,t-^C2fi''xi^Cifl''t (4.3) 

for n,m £ Z+. 



Next, let us discuss the explicit formula in the fourth quadrant. Naively, we use the initial 
conditions /o-i = c^t^'' and fio = cq to get the formula /„_,„ = fn.m\ci>^c3f-'-- However, this setting 
makes the discrete power function become a single- valued function on Z^. In order to allow 
fn,m to be multi-valued on Z', we introduce a discrete analogue of the Riemann surface by the 
following procedure. Prepare an infinite number of Z^ -planes, cut the positive part of the "real 
axis" of each Z^-plane and glue them in a similar way to the continuous case. The next step is to 
write the initial conditions (3.1) and (4.1) in polar form as 

f(l,nk/2) = Ckt"- (A: = 0,1,2,3), (4.4) 

where the first component, 1, denotes the absolute value of n + im and the second component, nk/2, 
is the argument. We must generalize the above initial conditions to those for arbitrary 6 Z so 
that we obtain the explicit expression of /„,„, for each quadrant of each Z^ -plane. Let us illustrate a 
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Figure 6: The discrete power function 
with y - 512 whose domain is 7?. 



Figure 7: The discrete power function 
with y - 512 whose domain is the discrete 
Riemann surface. 

3 

typical case. When -n < arg(n + im) < In, we solve the equations (2.7) and (2.3) under the initial 
conditions 

f{\,3nl2) = c,t"\ /(I, In) = c,t"\ (4.5) 

to obtain the formula 

f-n-m= fn,m\cQ^c^&,Ci^Cif-'^ («, m 6 Z+). (4.6) 

We present the discrete power function with 7 = 5/2 whose domain is 7? and the discrete Riemann 
surface in Figure 6 and 7, respectively. Note that the necessary and sufficient condition for the 
discrete power function to reduce to a single- valued function on 7? is (q = cm and) e'*""' = 1, 
which means that the exponent y is an integer. 

5 Associated circle pattern of Schramm type 

Agafonov and Bobenko have shown that the discrete power function for real y is an immersion and 
thus defines a circle pattern of Schramm type. We have generalized the discrete power function to 
complex y. It is natural to ask whether there are other cases where the discrete power function is 
associated with circle patterns. We have the following result: 

Theorem 5.1 The mapping f : ^ C ; (n, m) i-^ /„ ,„ satisfying (2.1) and (2.3) with the initial 
condition 

/o,o = 0, /i,o = l, hi=e'''^'\ yeC (5.1) 
is an immersion when Re y = 1. 

In this section, we give the proof of Theorem 5.1 along with the discussion in [J]. We also use the 
explicit formulae given in previous sections. 
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5.1 Circle pattern 

Setting 



y=l+iS, SeR, (5.2) 



we associate the discrete power function with circle patterns of Schramm type. The proof of 
Theorem 5.1 is then reduced to properties of the radii of those circles. 

Lemma 5.1 A discrete power function /„ „, defined by (2.1) and (2.3) with initial condition 

/o,o = 0, /i,o = l, /o,i = Cle'^'>/^ ci>0 (5.3) 
for arbitrary y e C\2Z has the equidistant property 

flnfi ~ fln-lfl = /2«+l,0 ~ flnfi, /o,2m ~ /o,2m-l = /o,2m+l ~ fo,2m, (5.4) 

for any n > l,m > 1. Moreover, if and only if Re y = 1, we have 

l/«+l,0 ~ fn,o\ = \fn,0 ~ /n-l,ok l/o,m+l ~ fo,m\ = l/o,m ~ /o,m-ll; (5.5) 

for any n > l,m > 1. 

Proof. By using the formulae in Theorem 3.1 (or Proposition 3.2), we have 

y t\ 2k + y 



flnfl — fln-lfi - fln+lfl — flnfi - FT 



J 0,2m - J 0,2m- i - J0,2m+l - J 0,2m - CyC ^ — | | _ 



(5.6) 



2n + y \ \ 2k- y 

k=\ 

niy y 2k + y 

which proves (5.4). We also have 

_ fr 2fc - 2 + y _ A 2fc - 2 + y 

j2n+2,0 - J2n+l,0 - — > J2n+l,0 - J2n,0 - — TT, • (.->■') 

Putting y = 1 + iS, we obtain 

_'f. 2k-l+i6 _f. 2k-l + i6 

j2n+2,0 - j2n+\,0 " | | 2k - \ - iS ' -'^n+l.O " J2n,0 - [ [ 2^ - \ - ig ' ^ ' 



k=\ k=\ 



which implies |/2„+2,o - /2n+i,ol = l/2«+i,o - f2n,o\ = 1- Using the first equation of (5.4), we see that 
the first equation of (5.5) follows. The second equation of (5.5) can be proved in a similar manner. 
Suppose that (5.5) holds, then from (5.7) we have 



/2n+2,0 ~ /2«+l,0 



/2;i+l,0 ~ /2n,0 

which leads us to Re y = 1 . □ 



2n + y 



2n + 2 - y 



1, (5.9) 
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Proposition 5.2 Let fn,m satisfy (2.1) and (2.3) in with initial condition (5.3). Then all the 
elementary quadrilaterals (fn,m, fn+\,m, fn+i,m+u fn,m+i) are of the kite form, namely, all edges at 
each vertex fn,m with n + m = \ (mod 2) are of the same length, 

\fn+\,m fn,m\ ~ \fn,m+l fn,m\ — \fn~l,m fn,m\ — \fn,m-l fi.ml- (5.10) 

Moreover, all angles between the neighbouring edges at the vertex fn,m with n + m = (mod 2) are 
equal to n/2. 

Proof. For three complex numbers Zi (i =1,2, 3), we introduce a notation 

[zi,Z2,Z3] = arg- — —. (5.11) 

We first consider the quadrilateral (/().o,/i.o,/i.b/o,i)- Notice that /qj = zcie^^ which implies that 
[/o,i,/o,o,/i,o] = f • Then it follows from (2.1) that [/i,o,/i,i,/o,i] = f , l/i,o - /o,ol = l/i,i - /i,ol and 
l/o,i -/o,ol = l/i,i -/o,il- We next consider the quadrilateral (/i,o,/2,o,/2,i,/i,i) where, from Lemma 
5.1, we have = |/2,o-/i,ol- We see from (2.1) that [/2,i,/i,i,/i,o] = [/i,o,/2,o,/2,i] = f and 

1/2,1 - = 1/2,1 - fifiV From Lemma 5.1 and [/i,o,/2,o,/2,i] = f , we see that [/2,i,/2,o,/3,o] = f • 
Then a similar argument can be applied to the quadrilateral (/2,o, /s.o, /3,i> /2,i) and so forth. In this 
manner. Proposition 5.2 is proved inductively. □ 

From Proposition 5.2, it follows that the circumscribed circles of the quadrilaterals 
(fn'i,m, fn.m-i, fn+i,m, fn,m+i) with n + m = 1 (mod 2) form a circle pattern of Schramm type [4, 16], 
namely, the circles of neighbouring quadrilaterals intersect orthogonally and the circles of half- 
neighbouring quadrilaterals with a common vertex are tangent (See Figure 8). Conversely, for 
a given circle pattern of Schramm type, it is possible to construct a discrete conformal mapping 
fn^m as follows. Let {Cn,,n}, (n,m) G V = {(n,m) e Z^\n + m = 1 (mod 2)} be a circle pattern of 
Schramm type on the complex plane. Define / : ^ C ; (n, m) i-^ /„ ,„ in the following manner: 

(a) If n + m = I (mod 2), then f„j„ is the center of C„_,„. 

(b) If n + m = 0(mod2), then f„^,„ := Cn-\,m n C„+i,,„ = C„,„,+i n C„,„,-i. 

By construction, it follows that all elementary quadrilaterals (fn,m, fn+i,m^ fii+i,m+u fn,m+i) are of 
the kite form whose angles between the edges with different lengths are |. Therefore, (2.1) is 
satisfied automatically. In what follows, the function f„ „„ defined by (a) and (b), is called a discrete 
conformal map corresponding to the circle pattern {C„ „,}. 

We now use the the radii of corresponding circle patterns to characterize the necessary and 
sufficient condition that the discrete power function is an immersion. 

Theorem 5.3 Let satisfying (2.1) and (2.3) with initial condition (5.3) be an immersion. Then 
Rnjn defined by 

Rn,m ~ \fn+i,m fn,m\ \fn,m+l fn,m\ \fn-l,m fn,m\ \fn,m-l fn,m\ (5-12) 

satisfies 

Rn+l.m Rn-\,m Rn,m+1 Rn,m-\ ^ /r- i i\ 

n h m = 0, (5.13) 

Rn+l,m "I" Rn—\,m Rn,m+\ Rn,m—l 
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T 

- /3.2 

T 

£ /3.1 



/3,0 



Figure 8: A schematic diagram of quadrilaterals. 

and 

[(m + l)i?„^m+l + 6Rn+\,m\Rn,m+\{Rn+l,m + Rn-l,m) + ,„+i — Rn+l,mRn-\,m) 

-^/7+l /77+2 2 ' 

[(m + l)R„+i,,n — 6R„^,n+l](Rn+l,m + Rn-\,m) ~ ^(Rn m+l ~ Rn+l,mRn-l,m) 

(5.14) 

for (n,m) e V. Conversely, let R : V ^ R+ satisfy (5.13) and (5.14). Then Rn^,„ defines an 
immersed circle pattern of Schramm type. The corresponding discrete conformal map fn,m is an 
immersion and satisfies (2.3). 

Proof. The proof of Theorem 5.3 occupies the remainder of this subsection. Suppose that the 
discrete power function fn^m is immersed. For n + m = 0(mod2), we may parametrize the edges 
around the vertex as 

fn+\,m fn,m ~ '"l^'^j fn,m+l fn,m — i^l^^j fn-l,m fn,m — ^3^^^ fn,m-\ fn,m — ir/\e^ , (5.15) 

where r, > (z = 1, 2, 3, 4) are the radii of the corresponding circles, since all the angles around 
fn^m are |. The constraint (2.3) reads 

yfn^m = e"" lln^^ + 2/m-^) . (5.16) 



■- h.2 



■'x/i.: 



£3-< 



y' /i. 



/2.1 



Lemma 5.4 For n + m = (modi) we have: 

f _ f - _ ^1 ~ '"^2 , _ , _ ip ri + ;>4 

/>!+l,m+l Jn+l,m — ^\ . ■> Jn+l.m Jn+l,m-l — ^ f\ . • 

ri + ir2 ri - ir^ 
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/i-ljn+1 f„.m+\ fn+\.m+\ f,i+2j„+l 



f,,-,., 



J 


L 

''1 ' 

ri /f 

g ^(-^ 


J - 

R, ■ 

r\ 


• E 

'■3 fn..„ 

> ri , 

n 


r 


E 

/„+2.,„ 



Figure 9: Parametrization of edges. 



Proof. The kite form of the quadrilateral {fn,m, fn+i,m, fn+im+i, fn,m+i) implies /„+i,„,+i - = 
-ifn+i,m - (See Figure 9). The first equation of (5.17) follows by noticing 

that tan[y;, „,+i,y;,+i „,,/„+! m+i] = ^. The second equation is derived by a similar consideration on 

the quadrilateral {fn,m-i, fn+\,m-\^ fn+l.m^ fn,m)- c 

Setting /„+2,m - /n+i,m = He'*^"^""^ and substituting (5.15)-(5.17) into (2.3) at the point (n + l,m), 
one arrives at 

2n^^ + 2im^^ + yn = 2(n + l)ri -+2mn —■ (5.18) 

ri + r2 + r4 I + e 2r\{r2 + r^) 

The real part of (5. 1 8) gives 

n h m = 0, (5.19) 

ri + r3 r2 + r4 

which coincides with (5.13). 

Now we parametrize the edges around the vertex fn+\,m+\ with n + m = (mod 2) as 

fn+2,m+\ — fn+l,m+l + ^1^^ , fn+l,m+2 — fn+l,m+l "t" Z^2^^ > 
fn,m+\ ~ fn+ l,m+l - /"l^ , i,m+i - ir\e . 



(5.20) 



From the first equation of (5.17) and noticing that all angles around the vertex fn+\,m+\ are f , we 
have the following relation between yS' and yS: 

e^P' = ie'l^HzlH. (5.21) 
n + ir2 

One can express fn+2,m+\ in two ways as (See Figure 9) 

/„+2,.+i = U2,m + iRye'^^''^ = Uu,n + ^e'^""^ + iRye'^^^'K (5.22) 

and 



y;,+2,,„+i =y;,+i,,.+i+i?ie'^' =/„+i,m-e'^i^^^^+i?i^^'^'. (5.23) 

ri + ir2 

The compatibility implies 



= ^iHnnzln. (5.24) 

Ri - iri ri + i>2 
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Then, from the imaginary part of (5.18), one obtains 



rirA - n 



m- 



+ 6ry ={n + 1)- 



riRi 



ri+Ri 



or solving (5.25) with respect to Ri = i?„+2,m+i 

5Rn.m+ 1 ]-^«+ 1 ,m(.Rn,m+\ 



(5.25) 



R 



n+2,m+l 



{{n + l)Rn,m+\ + Rn,n-i) - m{Rf,^^^^^^ - R 

n,m+\Rn,m-\) 

We may rewrite (2.3) at (n + 1, m + 1) in terms of r,- (/ = 1, 2, 3, 4) R2 as 



2n h 2im 1- yri 1 - 



(5.26) 



ji - ir2 



n + ir2 



(5.27) 



= i- 



n + ir2 

Eliminating 7 from (5.18) and (5.27), we get 

r2 



2(n+ 1) „ +li{m+ 1)- 



m- 



r2 + r4 ri + 
We then eliminate Ri using (5.25) to obtain 



n (m+1) —+in + l)- 



r2 + i?i 



= 0. 



(5.28) 



R7 



[(m + l)r2 + 6ri]r2(ri + r3) + nri(r^ - r^rj,) 
[(m + l)ri - 5r2](ri + r^) - n(rl - rir^) 



(5.29) 



which coincides with (5.14). This proves the first part of Theorem 5.3. 
To prove the second part, we use the following lemma. 

Lemma 5.5 LetR : V ^ K.+ satisfy (5.13) and (5.14). Then, R satisfies (5.26) and 

[(^ + l)(^Ll,;;i ~ Rn,m-lRn+2,m-\) + SRn+\,m(Rn.m-\ + Rn+2 ,m- 1 n,m+l + Rn,m-l) 

~f^iRl+l,m ~ Rn,m+\Rn,m~\){Rn.m-l + Rn+2,m-l) = 0. 



(5.30) 



Proof. Substituting (5.13) at (n, m) and at (n + 1, m+1) into (5.14) to eliminate i?„-i,„, and7?„+i,„,+2, 
we get (5.26). Substituting (5.14) at (n + l,m - 1) into (5.26), we get (5.30) under the condition 
5Rn+i,m{Rn+2,m-\ - Rn,m+i) + (n + M + l)(Rn+2,m-iRn,m+i + RI+i,„) ^ wWch cau bc Verified from the 
compatibility with (5 . 1 3 ) and (5.14). □ 

Eliminating 6 from (5.26) and (5.30), we get 



n+l,m 



-R„ 



n,m+lRn+2,m+\ 



)(Rn,m-l + Rn+2,m~l) + (Rn+\,m Rn,in—lRn+2,in- \)(Rn,m+i +Rn+2,m+l) - 0. (5.31) 



In [16], it was proven that, given R„ ,„ satisfying (5.31), the circle pattern with radii of the circles 
i?„ ,„ is immersed. Thus, the corresponding discrete conformal map yj,_,„ is an immersion. 

Let us finally show that the discrete conformal map satisfies (2.3). Putting m = in (5.13), 
we have Rn+i,o = R„-ifl. This means that |/„+i,o - fn.ol = \fn,o - fn~\fi\ for any n > 1. By using the 
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h-N.\ 



n D 























0' '7- 





si D 



/2N+1.0 

(i) < 20 < f 



flN.O 



flN +1.0 

2 (ii) f < 20 < TT 

Figure 10: Configuration of points in Lemma 5.6. 



ambiguity of translation and scaling of the circle pattern, one can set /o = 0, /i = 1 without loss 
of generality, and set 



Putting {n, m) = {IN, 0) in (5.26) we have 

5 = (2A^+ 1) 



( N \ 



(A^= 1,2,...). 



(5.32) 



V 7=1 / 



-^2^+1,0 ~ -^2^1-^2^+2,1 



(5.33) 



RlN+l.aiRlN,! + R2N+2,l) 

A geometric consideration leads us the following lemma: 
Lemma 5.6 We have 

5 = (2A^+ l)tan0;v+i (A^ = 0, 1, . . .). (5.34) 
Proof. First, we consider the case of < 20 < nil, see Figure 10 (i). From BC = AD we obtain 

(r2 +i?i)sin20' + (/z + 7?i) sin20 = 2ry, (5.35) 

which yields 



Since + r\> 0, we have h = 



(r2h + r^) [(r2 + Ri)h + (r2i?i - r])] = 0. 
- r2R, 



(5.36) 



r2+Ri 



and 



tan0 



- r2Ri 



(5.37) 



riir2+Ri) 

Thus we get (5.34) from (5.33). When n/2 <26 < n, the configuration of points is shown in Figure 
10 (ii). The equality BC = AD implies 

(r2+i?i)sin20' + (/z+7?i)sin(;r-20) = 2ri, (5.38) 

which gives the same result as the case of < 20 < n/2. Let us investigate the case of -7t/2 < 
20 < 0, see Figure 1 1 (i). Since DC = AB we have 



(h + Ri) cos 20 + (r2 + i?i) cos(;r - 20') = h + r2. 



(5.39) 
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f2N+2,\ 



(i) -f < 20 < (ii) -n<2e< -f 

Figure 11: Configuration of points in Lemma 5.6. 



which leads us to 

(r2h + r?) [(r2R, - r\)h - r\{r2 + Ri)] = 0. (5.40) 

Then we have h = —, and thus (5.37). Figure 11 (ii) illustrates the case of -n < 29 < 

riRi - r\ 
-nil. We see from AB = DC that 

{h + Ri)co%{n - \ie\) + h + r2 = (r2 + Ri)cos(n - 26'), (5.41) 

which also leads us to (5.40), and thus (5.37). Therefore we have proved Lemma 5.6. □ 

From (5.32), (5.34) and the initial condition /o,o = 0, /i,o = 1, we see by induction that the points 
fnfi satisfy 

r ^ (fn+lfi - fn,o)(fnfi - fn~l,o) ^CA^\ 

7 fnfi = 2n . (5.42) 

Jn+lfi — Jn- 1,0 

Similarly, we see that /o.,„ satisfy 

r ^ ifo.m+l ~ fo.m)ifo,m ~ fo,m-l) ,r- .^-^ 

yfojn = 2m — ■ . (5.43) 

/o,m+l ~ /o,m-l 

Thus it is possible to determine f^^^ in by using (2.1). Since (2.1) is compatible with (2.3), fn,m 
satisfies (2.3) simultaneously. This proves the second part of Theorem 5.3. □ 



5.2 Positivity of radii of circles 

Theorem 5.3 claims that if 7?„,„, satisfying (5.13) and (5.14) is positive, then the corresponding 
/„„, is an immersion. In order to establish Theorem 5.1, we have to prove the positivity of i?„„, 
determined by (5.13) and (5.14) with the initial condition 7? l o = 1, Rq,\ = ^(> 0). First, we show 
that positivity of for (n, m) e V is reduced to that of Rn,n+i for n 

Proposition 5.7 Let the solution i?„_„, of (5.13) for (n,m) 6 V and (5.14) for n = m e Z+ with 
initial data 

i?i,o=l, Ro,i=a>0) (5.44) 

be positive for m = n + \, namely, Rn,n+i > Ofor any n 6 Z+. Then Rn,m is positive everywhere in V 
and satisfies (5.14) for (n, m) e V. 
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Proof. Equation (5.14) for (n, m) = (0, 0) with initial data (5.44) determines R\2- We use (5.13) 
and (5.14) for n = m inductively to get Rn,n+i and Rn+i,n- As was mentioned before, we see that 
Rin+ifi = 1 and Ro^im+i = ^ for all n,m e Z+ by putting n = and m = 0, respectively, in (5.13). 
With these data one can determine i?„,„, in V by using (5.13). When n > m, we use (5.13) in the 
form of 

(n - m)R„^,„+i +(n + m)Rn,m-\ ^ 
{n + m)R„^,„+i +(n- m)R„^^.i 

For positive Rn-i.m^Rn,m+\ and Rn,m-\, we get i?„+i,„, > 0. When m> n, one can show in a similar 
way that i?„,„,+i > for given positive R„^,n-i,Rn+\^m and Rn-\,m by using (5.13). One can show by 
induction that we have (5.30) for m = n > I, and we get (5.14) for n = m + 2. Similarly, one can 
show by induction that we have (5.30) at {n + 2k, n) for n,k > 1, and (5.14) for {n + 2k, n). Thus 
we obtain (5.14) for n > m. One can show in a similar way that we have (5.14) for n < m by using 
(5.30) at (n, n + 2k) as an auxiliary relation. □ 

Due to Proposition 5.7, the discrete function with Re y = 1 is an immersion if and only if 
Rn,n+\ > for all n € Z+. We next reduce the positivity to the existence of unitary solution to a 
certain system of difference equations. 

Proposition 5.8 The map / : ^ C satisfying (2.1) and (2.3) with the initial condition fa o = 
0, /i,() = l,/o,i = i^i^ > 0) is an immersion if and only if the solution (x„,_y„) to the system of 
equations 



(Xn - 1) + (Xn + 1) — - = 0, 

J n + 2 n + I y 



(5.46) 



O - 1 -1-1 1,-1 ' ^0^0 ~ T 

2 l-x'y' l+^,/y„+i r-2 



with 



yo = (5.47) 



is of the form x„ = e^'"",yn = e^"^", where a„,(pn £ (0,n/2). 

Proof. Let be an immersion. Define a„, (fn e (0, n/2) through 

fn,n+2 ~ fn,n+l — ^ " {fn+l,n+l ~ fn,n-\-\\ fn+l,n+l ~ fn,n+l — ^ '^"i.fn,n ~ fn.n+l)- (5.48) 

Using Proposition 5.2, one obtains 

f , , - f , - R Ji2ifi„-7r/2+l3„) f , _ f , , - P J(J3„+2a„-i-7i/2) 

Jn+l,n+\ Jn,n+l ^^n,n+\'^ 5 jn,n+\ jn~\,n+\ ^^f!,«+l'' ; 

f ^ - f , - R J(2(f„-JT/2+2a„+/i„) f f _ n i(fi„+jT/2) 

Jn,n+2 Jn,n+l ^^Ji.n+l'- 5 Jn,n+l Jn,n ^^n,n+l'-- 5 



(5.49) 



where /3„ = arg (/„+i_„ - /„,„). Figure 12 is a schematic diagram illustrating the configuration of the 
relevant quadrilaterals. Now the constraint (2.3) for (n, n + 1) is equivalent to 

r/„..i = -2/i?„,.ie'^" ( ^_„,„ + ^-.,11,,) _i ) • (5.50) 

Putting these expressions into the equality 

fn+l,n+2 — fn,n+\ = s'^^"^" '^^'^^^"\Rn,n+\ + iRn+l,n+2), (5.51) 
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fn-l,n+2 fn,n+2 fn+\,n+2 



L 

9 Rn,n+\ - 

fn,n+i 


J 

a„ / 

"X 01,1 


3 ^ 7 

n . r 


E 

^\ Vn Rn,n+1 
1 



Vn+l,n+2 



fn-l,n fn,n fn+\,n 

Figure 12: A schematic diagram of quadrilaterals for Proposition 5. 

together with i?„+i,,i+2 = Rn,n+\ tan cr,, and = e'(/3"-'^/2+2y«)^ one obtains 

« + 1 n + 1 



2i sin Qf„ 
+2 cos a„ 



n n+1 



+ 



]^ -I- g2!(i/'„-ff„-i) ' g-2m„ _ g2(V„ 

On the other hand, equation (5. 14) for m = n is reduced to 



SRn,n+l =(n+ 1)- 



• n- 



By a similar geometric consideration to the proof of Lemma 5.6, this implies 

S = -in + 1) cot(Qr„ + ipn) - n tan(a„_i - 
which can be transformed to 



7 



n + 1 



(5.52) 



(5.53) 



(5.54) 



(5.55) 



By using (5.55), we see that (5.52) yields the first equation of (5.46) with jc„ = e^'"" and y„ = e^"^". 
The second equations of (5.46) come from (5.55). This proves the necessity part. 

Now let us suppose that there is a solution (x„,j„) = {e^^"" , e^"^") of (5.46) with an,(Pn 6 (0,7r/2). 
This solution together with (5.48) and (2.1) determines a sequence of orthogonal circles with their 
centers on and thus the points {fn,n+\, fn±\,n^i, fn,n, fn^n+i)- Now (2.1) determines on 

Z^. Since q'„,^„ € (0,;t/2), the inner parts of the quadrilaterals (fn,n+i, fn+\,n+u fn+i,n+2, fn,n+2) and 
of the quadrilaterals {fn,n, fn+i,n, fn+\.n+i, fn,n+\) are disjoint, which means that we have positive 
solution Rn,n+i and Rn+i,n of (5.13) and (5.14). Given i?„+i,„ and (5.13) determines 7?„ ,„ for 

all {n, m) e V. Due to Proposition 5.7, 7?„ „, is positive, and satisfies (5.13) and (5.14). Theorem 5.3 
implies that the discrete conformal map gn,m corresponding to the circle pattern {C,j m} determined 
by Rn,m is an immersion and satisfies (2.3). Since gn,,i = fn,n and gn,n±\ = fn,n±\, equation (2.1) 
implies g„jn = fn,m- This proves Proposition 5.8. □ 
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Note that although (5.46) is a system of equations, a solution (jc„,y„) of (5.46) is determined by its 
initial value 3^0. 

The system of equation (5.46) can be written in the following recurrent form: 

l+xl- IXnJn _ (2n + 2 + y)Xn + yyn+\ ,^ 

or 

y„+i = ^{Xn,yn), Xn+i ='i'{n,Xn,y„+i), (5.57) 

x-y +xy-'-2 l (2n + 3 + /g).yy-'+(l+''^) 

xy + X - 2 x (2n + 3 - zo).x; + (1 - zo) 

and xq = _ ^ J_ j|_ |g g^gy |<l)(.x:,y)| = \^(n,x,y)\ = 1 when |x| = lyl = 1 and that 

l-i6yo 

\xq\ = 1 when \yo\ = 1, which implies that this system possesses unitary solutions. Moreover, we 
have the following theorem as for the arguments of the unitary solutions: 

Theorem 5.9 There exists a unitary solution (Xn,yn) to the system of equation (5.46) with Xn,yn e 
Ai\{+1}, where 

Ai = {e2'^|;8e[0,;r/2]}. (5.59) 

Proof We first investigate the properties of the function 0(x, y) and *F(n, x, y') restricted to the 
torus = 5^ x5> = {{x,y)\x,y e C, |jc| = \y\ = 1}. 

Property 1. The function (^{x,y) is continuous on Aj x Ai\{(+1, +1)}. The function ^'{n,x,y) is 
continuous on Aj x Aj for any n e Z+. (Continuity on the boundary of Aj x Aj is understood to be 
one-sided.) 

The points of discontinuity must satisfy 

x^y + y-2x = 0, (2n + 3 - i5)y + (1 - i6)x = 0. (5.60) 

The first identity holds only for ix,y) = (+1, +1). The second never holds for unitary x,y. 

Property 2. For (x,y) e Aj x Ai\{(+1, +1)}, we have 0(x,y) 6 Aj. For (x,y) e Aj x Aj, we have 
^(n, x,y) G Ai U An U Ajv, where An := {e^'^ \/3 e (n/2, n]} and Aiv := {e^'^ \/3 e [-n/2, 0)}. 

Property 2 is verified as follows: using the transformation 

Un = tana„, v„ = tan^„, (5.61) 

where x„ = e^'""and yn = e^"^", we see that the first equation of (5.46) takes the form 

v„+i = u~\n- (5.62) 

It is obvious that m~^v 6 [0, +oo] when (u, v) 6 [0, +oo] x [0, +oo]\{(0, 0), (+oo, +oo)}. The second 
equation of (5.46) can be expressed as 

n 2,Y. (2n + 3 + i6)xny;,li + (I + iS) 
an^i = ojn -an + -, fi''"" = _^ — — . (5.63) 

2 {2n + ?,-i5)x„^yn+\+{l-i6) 
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By using the variables u,,, Vn+i, we have 

tan(a>„ - a„) = F(n, Un, v„+i), (5.64) 

where 

, [(n+\)-(n + 2y]u-(2n + 3)v + S(\+uv) ^. _ 

r{n,u,v) = . (5. 65) 

^ (n + 2) - (n + l)v^ + (In + 3)uv + Sv(l + uv) 

Lemma 5.10 It holds that cOn - or,, + | e [-7r/2, n] for a„, e [0, n/2]. 

Proof. Let us investigate the function F{n, u, v) for w, v e [0, +oo]. It is easy to see that 

dF(n,u,v) dF{n,u,v) 

> 0, < (5.66) 

ou ov 

on [0, +oo]2 except for the points satisfying (n + 2) - (n + l)v^ + (2n + 3)uv + 6v(l + uv) = 0. 
Consider the values of F(n, u, v) on the boundary of [0, +00]^. It is easy to see that 

F(n, 0, +00) = 0, F(n,0,0) = — F(n, +00, +00) = F(n, +00, 0) = +00. (5.67) 

n + 2 6 

dF{n,0,v) 6+ JS^+A{n + l){n + 2) 

We find that < except for the point v = v+ := , and that 

dv 2{n + 1) 

{n + l)u + 5 

F(n, u, 0) = IS monotone increasing. Note that 

n + 2 

(n + 2)u -(n + 2y + 6v + n + I 

F{n, u, +00) = — — — — , F(n, +00, v) = — — . (5.68) 

n + I — ou v(ov + 2n + 5) 

dF(n,u, +00) . n+l , , dF(n,+oo,v) 

When d > 0, we see that > except for the point u = and that < 

du 6 dv 

dF{n,u,+oo) JF(n, +oo,v) 

except for the point v = 0. When 5 < 0, we see that > and that < 

du dv 

2n + 3 

except for the points v = 0, := . The singular points of F(n, u, v) can be expressed by 

6 

{n + ly -6v-(n + 2) 

(u,v) = (G(v), v), where G(v) = , if v ?^ 0, v*. Note that (m, v) e [0, +00Y. 

v(ov + 2n + 3) 

Then we find that the singular points lie in v e [v+, +00], u e [0, {n + I)/ 6] (when 6 > 0) or in 
V G [v+,Vh.],m e [0, +00] (when 6 < 0), and that G(v) is monotone increasing (See Figure 13). 
Therefore we see that a*,, - a,, 6 [-n, n/2] when a,,, (p„+i e [0, n/2]. □ 

The final equation in (5.46) leads us to 

ao = -ifo + ^ + tan 6* = (5. (5.69) 

It is easy to see that -<pq + 9 + n/2 6 [-n/2,n] when (po e {0,n/2). Therefore property 2 is 
established. 

Now let us introduce 

Siiik) := {yo £Ai\xk£ An, xi 6 Aj (Z < k)], 
Sw(k) := {yo 6 Ai e Ajy, xi e Aj (/ < k)}. 
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+ 0O 



(0,0) (oo,0) (0,0) (oo,0) 

6>0 6<0 



Figure 13: Behaviour of m, v). 



wliere (x„,y„) is the solution of (5.46). From property 1, it follows that Sii(k) and Siv(^) are open 
sets in the induced topology of Aj. Denote 

Su = USuik), 5iv = U5iv(A:), (5.71) 



which are also open. These sets are nonempty since 5'ii(l) and 5iv(l) are nonempty (See the proof 
of Lemma 5.10). Note that 5ii(0) or 5iv(0) can be empty. Finally, introduce 

5i := {jo e Ai I x„ e Ai for all n e Z+}. (5.72) 

It is obvious that 5 1, 5 II and 5iv are mutually disjoint. Property 2 implies 

5iU5iiU5iv = Ai. (5.73) 



Since the connected set Ai cannot be covered by two open disjoint subsets 5ii and 5iv, we see that 
5 1 0. So there exists yo such that the solution (x„,);„) e Ai x Aj for any n e Z+. Suppose that 
an = 0, (fn 4^ hold at a certain n. Then we get = nil and a„+i = -nil from (5.56), or (5.62) 
and (5.63). Similarly, if = nil, ipn 4^ nil hold at a certain n, we get (^„+i = and then a^+i = n. 
It means that in both cases x„+i i Aj. Suppose that ^„ 4 0, ^„+i = hold at a certain n. Then we get 
a„ = nil and then Qr„+i = n from (5.56), or (5.62) and (5.63). Similarly, if 4 nil, (pn+i = ^11 
hold at a certain n, we get a,, = and then a„+i = -nil. It also means that in both cases Xn+i i Ai. 
Thus, it follows that a„ 4 0, nil and ip^ 4 0, nil for the solution {Xn,yn) 6 Ai x Aj for any n e Z+. 
□ 

We have shown that Si is not empty. In order to establish Theorem 5.1, let us finally show the 
uniqueness of the initial condition which gives rise to the solution (x„,j„) € Ai x Ai\{(+1, +1)}, 
namely, the circle pattern. Indeed, the initial condition is nothing but that for the discrete power 
function. Take a solution {x„,yn) such that y^ e Si and consider the corresponding circle pattern. 
Let Rk be radii of circles with centers at f2k,\, i.e., Ru '■= R2k,\- We have the following lemma. 



Lemma 5.11 An explicit formula for Ri_ is given by 

, . <-')'im-(f)j-'4m^(^)j 
'""(-i)'[m-mi-f[m-m]' 

where ^ = Rq. 
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Proof The radii Rk are defined by Rt = 1/2/1,0 - fikAl = |v2/t,ol|v2A.,i| (see (2.8)). From Proposition 
3.2, we have 

(r)k ipi-k+\,-r-k+l,-r + 2--\) 

V2k,0 - - — . 1^ ■, V2k,l 



where (p(a,b,c; -1) is given by (see (2.15)) 
(p{a,b,c;-l) = — — p{a,b,c;-l) 



(pi-k,-r-k+l,-r+l;-l) ' 



r(c) 



+ Cie 



Y{a-c+ l)Y{b - c + 1) 
r(2-c) 



F{a - c + l,b - c + l,2-c;-l), 



and 



^ = Cie-''^'\> 0). 
It is easy to see that \v2k,o\ = 1 under Rey = 1. Due to the formulae [1] 

T(a-b+ 1) 



F(a,b,a-b+\;-l) = 2-"n 



0^1/2 



Yi^a - b + l)Yi^a + 



F(a, b,a-b + 2;-\) = I'^'n^'^ib - l)"^r(a -b + 2) 

1 



X 



Y{\a)Y{\a - b + I) Y{\a + \)Y{\a - b + \) 



and the contiguity relation 



(c - a){c - b)tF(a,b,c+ \ ;t) + c[ia + b - 2c + l)t + c - 1] F{a,b,c; t) 

+ c(c- \){t- \)F{a,b,c- \;t) = 0, 



we get 



V2k.l = 



(-i)'[m+m]-'4m-(^)j' 

Then we can easily verify that V2kA satisfy the recurrence relation 

1 - iSkV2k,l 



V2(k+l),l 



V2k,l - l£k 



(5.75) 



(5.76) 



(5.77) 



(5.78) 



(5.79) 



(5.80) 



(5.81) 



(5.82) 



where Sk = 



2k + I 



On the other hand, (5.26) for m = is reduced to 



Rk+] — 



1 - SkRk 



Rk + Sk 

Then noticing R^ = ^, we see that Rk = \v2k,\\ = iv2k,\- □ 



Proposition 5.12 The set Si consists of only one element —3^ 



^-n6l2 ^ • 



■612 



(5.83) 
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flk.O flk+l.O flk+lfi 

Figure 14: Configuration of points around f2kfi for large k. 



Proof By using the formulae 



T{z)T{\-z) = 



n 



sin:7rz 



and the asymptotic formula as z ^ oo 



Y{az + ^) ~ {2nf'^e-"\azr^^-'^ (I arg z| < ;r , a > 0), 



we get 



_ (1 - + ^(1 + £-^-5/2) ^ (J ^ ^-nsn^^ - f(l - e-""^!^) 



(5.84) 



(5.85) 



(5.86) 



The relation (5.34) implies lim 6^ = 0. Figure 14 illustrates the configuration around fjko for 

sufficiently large k. Note that the three points f2k,o, fik+ifi, fik+i,o are asymptotically coUinear. 
We see from (5.86) that Rk should behave as lim i?^ = 1 irrespective of the parity of k, since the 

four points f2k+\,\, fik+iA, fik+i,!^ f2k+\,2 form a quadrilateral. Then we find that ^ = e'"^^^. This 
completes the proof of Proposition 5.12. □ 

Therefore f„^,n satisfying (2.1) and (2.3) is an immersion if and only if the initial condition is 
given by (5.1). This completes the proof of Theorem 5.1. 



6 Concluding remarks 

The discrete logarithmic function and cases where y e 2Z were excluded from the considerations 
in the previous sections. From the viewpoint of the theory of hypergeometric functions, these cases 
lead to integer differences in the characteristic exponents. Thus we need a different treatment for 
the precise description of these cases. However, they may be obtained by some limiting procedures 
in principle. In fact, Agafonov has examined the case where y = 2 and y = by using a limiting 
procedure [2, 3], the former is the discrete power function and latter is the discrete logarithmic 
function. In general, one may obtain a description of these cases by introducing the functions f„^m 
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and fn,,n as 



and 



fn,m • 



l (-r+l),. ^ 

lim. — — — /„,,n, 
r-'-J (r + l)y 



lim 

r->0 



forr = 2je2Z>o 
for 7 = -2j 6 2Z<o 

- 1 



(6.1) 



(6.2) 



respectively. The function /„ „, might coincide with the counterpart defined in section 6 of [4]. 

Moreover, it has been shown that the discrete power function and logarithmic function associ- 
ated with hexagonal patterns are also described by some discrete Painleve equations [5]. It may be 
an interesting problem to construct the explicit formula for them. 



Acknowledgements 

The authors would like to express our sincere thanks to Professor Masaaki Yoshida for valuable 
suggestions and discussions. This work was partially supported by JSPS Grant-in-Aid for Scien- 
tific Research No. 21740126, 21656027 and 23340037, and by the Global COE Program Education 
and Research Hub for Mathematics -for-Industry from the Ministry of Education, Culture, Sports, 
Science and Technology, Japan. 



A Backlund transformations of the sixth Painleve equation 

As a preparation, we give a brief review of the Backlund transformations and some of the bilinear 
equations for the r functions [ ]. It is well-known that Pyi (2.12) is equivalent to the Hamilton 
system 

, dH , dH , d 
op oq at 

whose Hamiltonian is given by 

H = khfAfi - [0^4/0/3 + «3/o/4 + (ffo - l)/3/4]/2 + ff2(ai + a2)/o. (A.2) 
Here ft and a, are defined by 

fo = q-t, h = q-l, f4 = q, f2=P, (A.3) 

and 

ao = 0, ai = Kao, as = /ci, 04 = kq (A.4) 
with ao + ai + 2a2 + + 04 = I. The Backlund transformations of Pyi are described by 

Siiaj) = aj - Uipi (/, J = 0, 1 , 2, 3, 4), (A.5) 
siifd = fi + ^, sif2) = f2-^ (i = 0, 3, 4), (A.6) 

J2 Ji 
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S5: ao<r^ au ^ a4, — — - — , /4 t—, 

'i' - JO 

r hihfi + 0^2) , t 
56 : ao ^3, ai a4, /2 , /4 1-^ — , (A.7) 

sj : ao<r-> a4, <r-> a^, f2 ^ , /a ^ — , 

f - ^ h 

where A = (aij)^j^f^ is the Cartan matrix of type \ Then the group of birational transformations 

{so, ■ ■ ■ , 57) generate the extended affine Weyl group IVCD^ '). In fact, these generators satisfy the 
fundamental relations 

^2 = 1 (/ = 0,...,7), SiS2Si = S2SiS2 (/ = 0,1,3,4), (A.8) 

and 



■^5>^{0,l,2,3,4| - 1,0,2,4,31 ■5'5, ■5'6>5'{0, 1,2,3,4] - ■5'{3,4,2,0,1|>5'6> ■5'7>5'{0,1,2,3,4| - ■5'{4,3,2,1,01>5'7, 



(A.9) 



We add a correction term to the Hamiltonian H as follows, 

Hq = H + + 4aiQ'2 + 4-al - {a^ + a4)^] + ^ [(ai + 0-4)^ + (as + 04)^ + 40204^ . (A. 10) 

This modification gives a simpler behavior of the Hamiltonian with respect to the Backlund trans- 
formations. From the corrected Hamiltonian, we introduce a family of Hamiltonians /z, (i = 
0,1,2,3,4) as 

t t-\ 1 
ho = Ho + -, hi = SsiHo) —, = S(,{Ho) + h4 = Sj(Ho), /zj = /ii + 5i(/zi). (A.l 1) 

Next, we also introduce t functions t, (z = 0, 1,2,3,4) by hi = (logr,)'. Imposing the condition 
that the action of the s^s on the t functions also commute with the derivation ', one can lift the 
Backlund transformations to the r functions. The action of ^^(D^ ^) is given by 

*o(To) = fo — , 5i(Ti) = — , 52(^2) = Al2lill!i^ ^3(^-3) = /3 — , S4{T4) = /4 — , (A.12) 
To Ti y/t T2 T3 T4 



(A.13) 



and 

55 : To 1-^ \t{t - l)]^Ti, Ti {t{t - l)]"5To, 

T3 r3(r - 1)?T4, T4 - 1)^?T3, T2 ^ [tit - l)]"5/oT2, 

i _i _i i _i 

56 : To 1-^ ?V'*r3, T3 1-^ -zY ^tq, ti i-^ ? ■^T4, T4 1-^ r^Tj^ T2 1-^ ? V4T2, (A. 14) 

57 : To (-l)"5(r - 1)^T4, T4 (-1)5(/' - l)~5To, 

Ti ^ (-l)ia - l)-ir3, T3 ^ (-l)-^a - l)'Ti, (A.15) 
T2 ^ -i(t- 1)~V3T2- 

We note that some of the fundamental relations are modified 

*,-*2(T2) = -525,(T2) (z = 5, 6, 7), (A.16) 
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and 

■^5 ■^6^(0, 1,2,3,4} = {i,-i,-i,-i,i}S6S5T{0,lX^A}^ 

^5^7T|0,1,2,3,4| = {i, -i, -1, -«>7^5T{0,1,2,3,4|, (A. 17) 

■^6'^7'^{0,1,2,3,4) = {-h-i,-^,i,i}S7S(,T[oi2,3A]- 

Let us introduce the translation operators 



(A. 18) 



(A.19) 



(A.20) 



Ti3 = 5152*0*4*2*1*7, 7^40 = *4*2*1*3*2*4*7, 
7^34 = *3*2*0*1*2*3*5, 7^14 = *1 *4*2*0*3*2*6, 

whose action on the parameters a = (ao, Ci, Q'2, Ci3, 0:4) is given by 

fi3((?) = (? + (0,1, 0,-1,0), 
f4o(a) = a + (-l,0,0,0,l), 
f34(a) = a + (0, 0,0, 1,-1), 
ri4((?) = (? + (0,1,-1,0,1). 

We denote Ti, i ,„ n, = T'^^T'"^T^^qT\^{to) {k, I, m, n' 6 Z). By using this notation, we have 

To,0,0,0 = To, T_i,_i,_i,o = [t{t- l)]^Ti, 

T0,-1,-1,0 = {-IT'^r^T^, To,-i,o,0 = (-1)~^(? - 1)^T4, 

T-1,-2,-1,1 = (-1)"5*o(To), To,-1,o,1 = {-lyhtit - 1)]?5i(Ti), 

T-1,-1,0,1 = -/n*3(T3), T_i_i,_ij ={t - 1)^54(^4), 

for instance. When the parameters a take the values 

{aQ,ai,a2,aj,,aa) = {-b,a + n' , -n',c - a,b - c + \ + n'), (A. 21) 

the function r^./,,,,,' relates to the hypergeometric r function r^/"' introduced in Proposition 2.6 
by[13] 

_ kj,m -{a+b-c+2n'flA-{a-b-n'fl4+n'{bW)-n'{n'-\)l2(. i n(S+6-c+2j!')^/4+1/2 ( s. 

'k,l,m,n' — <^/t,/,m,n' ' „' ' V ~ ^ ) ■> {/\.Z.Z) 

where we denote a = a + k,b = b + I + 1 and c = c + m, and the constants 0Jk,i,m,n' = 0Jkj,m,n'(a, b, c) 
are determined by the recurrence relations 

(^k+lJ,mjOJk-lJ,m,i = i^(c - «)ii>^ /„,,-, 

cok,i+umjOJkj~i,m,i = -ib(c - bWkx^j, = ^' ^) (A.23) 

(^k,i,m+i,iOJkXm-i,i = (c - a){c - b)a)^j^^^j 

and 

k,l,m,n' 



<^k,l,m,n'+l<^k,l,m,n'-l - ~^llmn' (A. 24) 



with initial conditions 



'^-1,-2,-1,1 = (-1)"^^"^^, ^^0,-2,-1,1 = b, 

0J-\-\~i,i = 1, ^0,-1,-1,1 = (-1)'^^'*, 

'^-1,0,0,1 = -(-l)^^'"^(c - a), ojQfifij = -i, 

_i,o,i = -i{c - a), ^^0,-1,0,1 = (-l)"^^'*, 
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(A.25) 



and 

<^-l -2,-1,0 = (-l)"^^"^^, 0^0,-2,-1,0 = -b, 

^-1,-1,-1,0 = 1> <^0,-l,-l,0 = (-1) ^ A 

_ f i\-3/4/' \ _ 1 

<i^-i,o,o,o - y-i^) yc - a), a>o,o,o,o - i> 
^i^-1,-1,0,0 = c-a, ii>o,-i,o,o = (-l)"^^"^. 

From the above formulation, one can obtain the bilinear equations for the r functions. For 
example, let us express the Backlund transformations S2(fi) = / + -7- = 0, 3, 4) in terms of the r 
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functions tj (j = 0, 1, 3, 4). We have by using (A. 12) 

a2r^T3T4 - 5i(Ti)525o(ro) + 5o(ro)525l(Tl) = 0, 

a2t~^ToT4 - 5i(Ti)52*3(T3) + *3(T3)*2*i(Ti) = 0, (A.27) 
a2t~^-ToTi - 5i(Ti)5254(T4) + 54(T4)52*i(Ti) = 0. 

Applying the affine Weyl group ^^(D^ ^) on these equations, we obtain 



4 

(ao + 0:2 + a4) r^T354(T4) - Si{Ti)S4S2So(Tq) + To5o>^4*2*i(Ti) = 0, 
(ao + 0:2 + a4) t^-TiT3 - T45452*o(To) + To5o52*4(T4) = 0, 

(ao + Ori + Qr2)r^T3T4 - Ti5i52*o(To) + To5o*2^i(Ti) = 0, 

(ffO + ffi + 0^2) nSi{Tx)T-i - 54(T4)5i52^o(To) + Tq^o* 1 *2*4(T4) = 0, 

(Q'2 + Q'3 + 014) r^ToTi - T4545253(T3) + T35352*4(T4) = 0, 

(^2 + 0^3 + «4) r554(T4)To - 5i(Ti)5452*3(T3) + T3535452*i(Ti) = 0, 



(A.28) 

(A.29) 
(A.30) 



and 



(A.31) 



a2r'-roT3 - 5i(Ti)5254(T4) + 54(T4)525i(Ti) = 0, 
(a\ + Q'4 + a2)f~'^TQT3 - T\SiS2S4{j4) + T4 54 52 ■^1(^1) = 0. 

For instance, the first equation in (A.28) can be obtained by applying sqS4 on the first one in (A.27). 
We also get the second equation in (A.28) by applying S()S4S(, on the second one in (A.27). Other 

■^n' 'T'm 'T'l 'T'k 
14^ 34^ 40 13 



equations can be derived in a similar manner. By applying the translation T"'.T^.T'T'^j to the 



bilinear relations (A.28) and noticing (A. 20), we get 

(OfO + a2 + a4-m)t 2T^,/-l,m-l,«'T*:-l,/-l,,„-l,«' + l 

+'''/l-,/-l,m,«' + l'!"/t-l,/-l,m-2,n' + TkJ,m,n'Tk-l,l-2,m-2,n' + \ = 0> 

(ao + a2 + <3'4 - m)Tk^ij-i^,n-l,n'Tk,l-hm-hn' 

+Tk,l~l,m,n'Tk~l,l~l,m-2,n' — Tk,l,m,n'Tk-\,l~2,m-2,n' = 0, 

and then (3.13) for the hypergeometric r functions. Similarly, we obtain for the hypergeometric t 
functions (3.14), (3.15) and (3.22) from (A.29), (A.30) and (A.31), respectively. The constraints 

fo = h-t, /3=/4-1, (A.33) 

yield 

To54*2*o(To) = 54(T4)*2*4(T4) - tT ^84828 ^{t y), 

To5i52*o(To) = r45l52*4(T4) - ?5i(Ti)52*i(Ti), (A.34) 
T354*2*3(T3) = *4(T4)*2*4(T4) -ri5452^l(Tl), 
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and 

T3^3(T3) - TqSq{Tq) = (t - 1)Ti5i(Ti), 
tTsSjiTs) - To^o(To) = (t- 1)T4^4(T4), 

from which we obtain (3.23) and (3.33), respectively. Due to (A.l 1) we have the relation 



ho-h = (t-l) 
Then we get the bilinear relations 



12/4 + ^(1 - ffs - 0-4) 



(A.35) 



(A.36) 



Dto • T3 = r-S4(r4)S2Si(Ti) + -(1 - ^3 - a4)roT3, 
Dto-T3= t^T4S4S2Si(Ti) + ^(1 - tti, + a4)ToT3, 
DTo ■ T3 = P-TiSiS2S4(T4) + ^{Oq - ai)TQT3, 



(A.37) 



(dg df\ 

where D denotes Hirota's differential operator defined by D g ■ f = t\ — / - g — . By applying 

\dt dt 

the translation T"'^T!^^T'^qT'^^ to the first bilinear relation of (A.37), one gets 



1 



1 



D + - \a3 + a4 - k + I + n' 



Tk,l,m,n' ' TkJ-l,m-l,n' - -t^{t - 1) lTk-\,l-l,m-l,n' + iTk+ 



l,l,m,n'-\ 1 



(A.38) 

which is reduced to the first relation of (3.38). The second and third relations of (A.37) also yield 
their counterparts in (3.38). 
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